Neutron Brillouin scattering with pulsed spallation neutron source 
- spin-wave excitations from ferromagnetic powder samples - 
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Neutron Brillouin scattering (NBS) method was developed using a pulsed spallation neu- 
tron source, and the feasibility of NBS was demonstrated by observing ferromagnetic spin 
waves in Lao.sSro.2Mn03 and SrRuOa powders. Gapless spin-wave excitations were observed in 
Lao.8Sro.2Mn03, which were continuously extended to the lower scattering vector Q from previ- 
ous results using single crystals. The novel result is a well-defined quadratic Q dependence in the 
spin-wave dispersion curve with a large energy gap in SrRuC>3 indicating robust spin-orbit coupling. 
PACS numbers: 75.30. Ds, 75.40. Gb, 75.50.Gg 



In investigations of spin mediated materials physics, 
it is essential to experimentally determine microscopic 
interaction energies. For instance, an exchange interac- 
tion within the mean field approximation corresponds to 
the ordering temperature. This value is directly derived 
by the momentum dependence of spin waves in ordered 
states. Thus, spectroscopic investigations using thermal 
neutrons have been established as the most valuable and 
powerful experimental methods, although actual experi- 
ments have been performed mainly with single crystals, 
which require large volume and high quality samples from 
advanced spectrometers installed even at the most in- 
tense neutron sources. 

The current research trend in material sciences requires 
new materials of high performance as well as energy effec- 
tive functioning for device applications. To support this 
requirement for magnetic materials as one typical exam- 
ple, ferromagnetic semiconductors are a key material for 
future spintronics devices, and they have already been 
proven to be functional at low temperatures. However, 
now, research directs to the discovery of new materials of 
high performance at ambient temperatures. Obviously, 
these materials are required to possess a high ordering 
temperature (Curie temperature T c for ferromagnets). 
Therefore, we need to seek magnetic materials with a 
large spin- wave stiffness constant and a zero energy gap, 
which show that they are easily magnetized. Most of 
newly synthesized materials are chemically and struc- 
turally complicated such that large single crystals are 
difficult to be synthesized. Neutron spectroscopy experi- 
ments require the use of powder samples, and sometimes, 
low quality materials with impurity phases. 

The neutron Brillouin scattering (NBS) method is the 
most suitable for measuring the momentum dependence 
of spin-wave excitations from (000) in the forward direc- 
tion. Here, NBS is defined in comparison with the optical 
Brillouin Scattering. The principles of NBS are not new, 
although, until recently, such experiments have been lim- 



ited to a very small energy range, because intense neu- 
trons have been available in the thermal energy region. 
Sub-eV neutrons are now used high flux neutron sources. 
In this work, we demonstrate the first clear result of spin- 
wave excitations from two powder ferromagnets by us- 
ing a pulsed spallation neutron source: the colossal mag- 
netoresistance manganite Lao.8Sro.2Mn03 (T c = 316 K) 
[EH and SrRu0 3 (T c = 165 K) [|. Both perovskites 
essentially possess a cubic lattice structure with a tiny 
distortion. The former is a typical example of a next gen- 
eration magnetic device. The latter has recently become 
attractive because of an enhanced anomalous Hall effect 
due to robust spin-orbit coupling of the Ru 4d orbital 
4i]. Large single crystals are not always available and, 
in particular, SrRuOs has not been synthesized as yet. 
However, powder samples are readily available. Scat- 
tering intensities from powder samples rapidly decrease 
as the scattering vector Q increases because of the pow- 
der average of the dynamical structure factor. Scattering 
intensities from ferromagnets are strong only in the for- 
ward direction near (000). Owing to the kinematic con- 
straints of neutron spectroscopy, incident energy {Ej) in 
the sub-eV region is necessary to measure scattering in 
a meV transferred-energy (w) range, and scattered neu- 
trons need to be detected at very small scattering angles 
(4>) of 1° or less, with a high energy resolution of Aui/Ei. 

Efforts for realizing such experimental conditions in 
order to apply the NBS method have been made at the 
High Resolution Chopper Spectrometer (HRC) installed 
at the Material and Life Science Facility (MLF) in the 
Japan Proton Accelerator Research Complex (J-PARC) 
Similarly, Pharos, at the pulsed spallation neutron 
source at the Los Alamos Neutron Science Center (LAN- 
SCE) in the Los Alamos National Laboratory [7| and 
BRISP at the steady state neutron source in the Institut 
Laue-Langevin (ILL) Q have been applied NBS. On the 
HRC, polychromatic neutron beams extracted from a 25 
Hz pulsed neutron source were monochromatized by a 
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FIG. 1. Scan loci with Ei = 102 meV at the centers of PSDs, 
which are located at the scattering angles from 0.5° to 2.8° 
(from left to right) in steps of 0.15° , used in the present exper- 
iment on the HRC. The (Q,u>) region above the dashed line 
(higher transfer energies) can never be accessed by a spec- 
trometer with the lowest scattering angle of (f> — 5°. 



Fermi chopper. The energy of the monochromatic beams 
can be adjusted to the type of measurement. The neu- 
trons scattered by the sample were then detected, and 
the energy transfer (ui) was determined from the time- 
of-flight of the detected neutrons. For NBS, scattered 
neutrons were detected by a two-dimensional array of 
position sensitive detectors (PSDs) at a small <j>. In the 
array, linear PSDs were mounted with their linear axes 
vertical, and the centers of the PSDs were located hori- 
zontally at <f> between 0.5° and 2.8° in constant steps of 
approximately 0.15° |6J]. The horizontal divergence of the 
neutron beams was controlled by a collimator mounted 
immediately upstream of the sample. The collimator was 
composed of slits of vertical cadmium sheets with a colli- 
mation of 0.3°. The measurement was performed at the 
temperatures T = 6 K and 254 K for Lao.sSro^MnOa 
and at T = 7 K for SrRuOa. Since an estimation of back- 
ground contribution was important, an empty scan under 
exactly the same geometry as the sample scan, with only 
the sample crystal removed, was also performed. 

The scan loci of the centers of the PSDs in the array 
on the HRC with Ei = 102 meV selected in the present 
experiment on the (Q, uj) space are presented in Fig.l. A 
detector element located at <f> scans points in the (Q,u>) 
space as described by 
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where H and m are Plank's constant divided by 27r and 
the neutron mass, respectively. It should be noted that 
low angle detectors are essential to access the present 
(Q,ui) space. In fact, the region above the dashed line 
in Fig. 1, which is ui = (h 2 Q 2 /2m)/ sin</>, or the en- 
velope of scan loci for cf> — 5° in eq. (1) with respect 



to Ei, can never be accessed by a spectrometer having 
the lowest scattering angle of = 5°, for instance. The 
high energy-resolution Au; /Ei, which was a few percent 
or less, was achieved with short-pulse neutrons extracted 
from a decoupled moderator (a neutron source for high 
resolution experiments) and delivered to the HRC. The Q 
resolution can be derived as follows from the derivatives 
of U) and <fi in eq. (1): 
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where h 2 k 2 /2m — Ei. The angular resolution for a small 
<j) is given by 
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where the symbols w and h indicate the width and 
height, respectively, m, s, and d denote the moderator 
(neutron source), sample and detector, respectively; and 
L' _1 = L^ 1 + L^ 1 , where L\ = 15 m and Li — 5.2 m 
are the distances from neutron source to the sample, and 
from the sample to the detector, respectively. Therefore, 
we can calculate AQ with the following parameters: w m 
= 100 mm, h m = aw m , w s = h s — 30 mm, and Wd = 12.7 
mm. The HRC faces the moderator area of 100 mm x 100 
mm, and an intensity gain of a factor of 3 was obtained 
for Ei — 102 meV by a guide tube. The collimator cut 
the intensity gain in the horizontal direction, because the 
collimation was almost identical to the beam divergence 
between the moderator and the sample without the guide 
tube. Therefore, a — 3 1 / 2 was used. For the plot of hd = 
92 mm in Fig. 2, the angular resolution was calculated to 
be Acj) = 17 mrad. In the present experiment, since the 
first term in eq. (2) was much smaller than the second 
term because of Alo determined below, the Q resolution 
could be approximated as AQ = hA(f> = 0.12 A -1 . Fig- 
ure 2 presents the energy spectra observed at each fixed 4> 
on the scan locus determined in eq. (1) after subtracting 
the empty scan data, where the neutron intensities were 
summed up over a detector length of 92 mm (±46 mm 
from the center position). Each spectrum shown in Fig. 
2 was fitted to the following spectral function: 
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for lo > meV by parameterizing A, B, ujq, T, and Aw, 
where n(w) + l = [1 — exp(— u/ksT)} -1 with ks being the 
Boltzmann constant. For uj < meV, the temperature 
factor n(oj) + 1 in eq. (4) was replaced by the product of 
the temperature factor and the detailed balance factor, 
n(uj) exp(w/A;BT). The magnetic form factor is denoted 
by F(Q) [10]. The first and second terms in eq. (4) cor- 
respond to a spin wave component and elastic scattering 
with the energy resolution width Aw, respectively. The 
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FIG. 3. Energy spectra for SrRu0 3 at T = 7 K obtained by 
a constant Q cut, which are integrated intensities within the 
range Q± 0.024 A -1 of the intensity map in the (Q,ui) space. 
The dashed lines are fitted curves to the observed spectra and 
the solid lines are spin wave components (see text). 



FIG. 2. Energy spectra for Lao.8Sr .2Mn0 3 at T = 6 and 245 
K observed on the scan loci at 4>- The dashed lines are fitted 
curves to the observed spectra, and the solid lines are spin 
wave components (see text). 

spin wave component is thus presented by the parame- 
ters A, uj and T. The peak intensity parameter A was 
equivalent for both T = 6 and 245 K: A = 0.28 ± 0.07 
and 0.26 ± 0.01 in the units in Fig. 2 for (j> = 0.64°, re- 
spectively, for instance, and therefore, the intensity dif- 
ference between T = 6 and 245 K can be explained by 
the temperature factor. The elastic component consists 
of a Bragg tail around (000) and incoherent clastic scat- 
tering. In fact, the width of the plot of B as a function 
of Q, which was converted for an appropriate <j> and with 
uj = meV from cq. (1), was obtained at 0.11 A -1 at 
the full width at half maximum by fitting to a Gaussian 
function centered at Q = A -1 plus a constant. The 
observed width was consistent with the above-calculated 
value of AQ. The value of B was constant for Q > 0.3 
A -1 , where the elastic component almost consists of in- 
coherent clastic scattering. The energy resolution was 
obtained to be Acv = 2.0 meV, i.e., Au/E, = 2%. 

Because every position in any PSD can be converted 
into an appropriate <fi, an intensity map in the (Q,co) 
space can be obtained by averaging the intensities on a 
ring, which require to all have the same Q, by using all 
the areas of the PSD array. Figure 3 presents the energy 
spectra spanned on the (Q, u>) space for SrRuOa at T — 
7 K. It corresponds to constant Q cuts, which were inte- 
grated intensities within the range Q± 0.024 A" 1 . Each 
spectrum was well fitted to eq. (4) for u> > meV and 
the equation in which n(u>) + 1 in eq. (4) was replaced by 
n(uj)cxp(uj/kBT) for uj < meV. The peak positions of 
the spin wave component were obtained from the param- 
eters wo and T. As shown in Fig. 3, the peak intensity 
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FIG. 4. Spin wave dispersion curves of SrRu03 and 
Lao.8Sro.2Mn03. The solid lines are fitted dispersion curves. 
The dotted lines are the boundaries of the present experi- 
ment (4> — 0.5° and 2.8°), and the dashed line is the upper 
boundary accessible with a spectrometer having the lowest 
scattering angle of <f> = 5° (see Fig.l). 

of the spin wave was nearly constant in Q, which was 
consistent with the powder averaged dynamical structure 
factor for ferromagnetic spin-waves, because the disper- 
sion relation was approximated with a quadratic form 
as described below. Similarly, the peak positions of the 
spin wave component were obtained from the constant Q 
cuts for Lao.8Sro.2MnOa; they were also obtained from 
the energy spectra on the scan loci at corresponding 
for SrRuC>3. The peak positions of the spin waves ob- 
tained from two independent analyses fall onto the single 
dispersion curve for each scan, as shown in Fig. 4. 

The spin wave dispersion curves in Lao.sSro^MnOs ob- 
served at T = 6 and 245 K were well fitted to a quadratic 
form lo = DQ 2 with D = 126 ± 15 and 91 ± 5 meVA 2 , 
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respectively. In previous works measuring spin waves us- 
ing single crystals, a dispersion curve without a spin gap 
with D = 131 meVA 2 (T = 14 K) and 89 meVA 2 (T 
= 250 K) for La . 8 Sr . 2 MnO 3 was reported (HQ. At 
T = 14 K and Q < 0.3 A -1 from an appropriate re- 
ciprocal point, the dispersion curve along [001] was well 
fitted to u> = DQ 2 with D = 131 meVA 2 and was iden- 
tical to that along [111] [2]. Therefore, the spin wave 
dispersion curves observed in the present powder sample 
were in good agreement with those in previous studies 
and confirmed the isotropic dispersion relation for small 
Q. Therefore, the NBS experiment proved effective for 
using the HRC in the research and development of high 
performance ferromagnetic materials. 

The spin wave dispersion curve of SrRu03 observed at 
T = 7 K was well fitted to a quadratic form with an en- 
ergy gap: u) = Eq + DQ 2 with E = 2.3 ± 0.6 meV and 
D = 47 ±8 meVA 2 . As mentioned above, the crystal lat- 
tices for SrRuOa and Lao.gSro.2Mn.O3 were assumed to 
have the same simple cubic structure so that the 6 neigh- 
boring magnetic ions were distributed along the three 
principal axes. In the mean field approximation, the 
Currie temperature is expressed as T c = 2zS(S + 1) J/3 
with z — 6 being the coordination number. The stiff- 
ness constant of spin waves is given by the linear spin 
wave theory of the ferromagnetic Heisenberg model to 
be D = 2zSJa 2 , with lattice constant a being the dis- 
tance between neighboring magnetic atoms. The D value 
obtained in SrRuC>3 corresponds to T c = 151 ± 26 K in 
the mean field approximation, which was in good agree- 
ment with the actual Curie temperature (T c = 165 K). 
Such a scaling between D and T c has been commonly ob- 
served in metal ferromagnets fiol ] . A ferromagnet having 
a high symmetry of crystalline structure usually shows 
gapless spin- wave excitations due to negligible single ion 
anisotropy, and, in fact, spin waves in Lao.sSro^MnOs 
exhibit a gapless dispersion relation. We clearly ob- 
served a large energy gap in the spin wave excitations in 
SrRuC>3, as described herein. This characteristic feature 
originates from an induced spin anisotropy, possibly from 
robust spin-orbit coupling acting on the Ru 4d orbital [1] . 
The large anomalous Hall effect, changing the sign in the 
middle between T c and the base temperature in the ferro- 
magnetic phase of SrRuC>3 was theoretically interpreted 
as a strong spin-orbit interaction This theory also 
predicted a magnetic monopole scenario, which must be 
proven by future experiments. Spin wave excitations in 
the present low Q region reveal a robust spin-orbit cou- 
pling effect in this crystal for the first time. 

In conclusion, we developed an NBS method using 
a pulsed spallation neutron source and demonstrated 
its feasibility by observing ferromagnetic spin waves in 



SrRuC-3 and Lao.sSr .2Mn03 powders. It should be em- 
phasized again here that the short pulse and high inten- 
sity of the neutron source in a sub-eV region, as well as 
the performance of the HRC, with its high energy reso- 
lution and low background, made the NBS method pos- 
sible. It showed that powdered samples could be used 
for neutron spectroscopy. This demonstration also shows 
a well-defined quadratic Q dependence in the spin wave 
dispersion curve and a large energy gap in the SrRuOs 
ferromagnet, which demonstrates the importance of the 
spin-orbit coupling effect. 

This neutron scattering experiment was approved by 
the Neutron Scattering Program Advisory Committee of 
the Institute of Materials Structure Science, High En- 
ergy Accelerator Research Organization (Proposal Nos. 
2011S01 and 2012S01). S. I. and T. Y. were partially sup- 
ported by a Grant-in- Aid for Scientific Research (C), and 
Y. E. was supported by a Grant-in-Aid for Challenging 
Exploratory Research, both were from the Japanese Min- 
istry of Education, Culture, Sports, Science and Technol- 
ogy. This work was also partly supported by the Funding- 
Program for World-Leading Innovative R&D on Science 
and Technology (FIRST Program) from the Japan Soci- 
ety for the Promotion of Science. 



[1] Y. Endoh and K. Hirota, J. Phys. Soc. Jpn. 66, 2264 
(1997). 

[2] F. Moussa, M. Hennion, P. Kober-Lehouelleur, D. 

Reznik, S. Petit, H. Moudden, A. Ivanov, Y. M. 

Mukovskii, R. Privezentsev, and F. Albenque-Rullier, 

Phys. Rev. B 76, 064403 (2007). 
[3] for instance, L. Klein, J. S. Dodge, C. H. Ahn, G. J. 

Snyder, T. H. Geballe, M. R. Beasley, and A. Kapitulnik, 

Phys. Rev. Lett. 77, 2774 (1996). 
[4] Z. Fang, N. Nagaosa, K. S. Takahashi, A. Asamitsu, R. 

Mathieu, T. Ogasawara, H. Yamada, M. Kawasaki, Y. 

Tokura, and K. Terakura, Science 302, 92 (2003). 
[5] S. Itoh, T. Yokoo, S. Satoh, S. Yano, D. Kawana, J. 

Suzuki, and T. J. Sato, Nucl. Instr. Meth. Phys. Res. A 

631, 90 (2011). 

[6] S. Yano, S. Itoh, S. Satoh, T. Yokoo, D. Kawana, and T. 
J. Sato, Nucl. Instr. Meth. Phys. Res. A 654, 421 (2011). 

[7] R. A. Robinson, R. J. McQueeney, and T. M. Kelley, 
Physica B 241-243, 161 (1998). 

[8] D. Aisa, S. Aisa, E. Babucci, F. Barocchi, A. Cunsolo, 
A. De Francesco, F. Formisano, T. Gahl, E. Guarini, A. 
Laloni, H. Mutka, A. Orecchini, C. Petrillo, W.-C. Pil- 
grim, A. Piluso, F. Sacchetti, J.-B. Suck, and G. Venturi, 
J. Non-Cryst. Solids 352, 5130 (2006). 

[9] R. E. Watson and A. J. Freeman, Acta Crystallogr. 14, 
27 (1961). 

[10] Y. Endoh and P. Boni, J. Phys. Soc. Jpn. 75, 111002 
(2006). 



